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ABSTRACT 

The feasibility of adding a coolant 
gas to the nozzle flow of a colloid thrus- 
tor to improve the particle formation ef- 
ficiency is examined. The study is based 
upon a thermodynamic analysis of the noz- 
zle flow, where a simple physical model is 
used to describe the energy balance of the 
condensation process. The noncondensable 
coolant gas provides an additional heat 
sink for the latent heat given up during 
condensation. In cases applicable to col- 
loid thrustors, the analysis indicated 
maximum efficiencies of about 80 percent 
using very optimistic assumptions. The 
analysis also showed that there is an op- 
timum fractional mass flow of coolant gas 
that maximizes the particle formation ef- 
ficiency. Finally because of boundary 
layer and non-continuum flow effects in 
the nozzle, the minimum allowable nozzle 
exit pressure is about 0.1 torr, while at 
the same time, because of space charge ef- 
fects in the accelerator, the maximum al- 
lowable nozzle exit pressure is limited to 
about 0.01 torr. Since these conditions 
are mutually exclusive, unique flow 
schemes appear necessary if a mixed-flow 
condensation-type colloid thrustor is to 
be of practical value. 

INTRODUCTION 

This analysis is an extension of one 
presented in reference 1. In that analy- 
sis, in which only a condensable vapor was 
used in a condensation-type colloid thrus- 
tor, the particle formation efficiency was 
found to be less than 50$ even under the 
most favorable conditions. As shown in 
reference 2, colloid thrustors should op- 
erate with an overall efficiency of better 
than 50$ in the 2000-5000 second specific 
Impulse range to be useful in electric 
propulsion. In reference 1 it was found 
that the main factor preventing higher 
formation efficiencies was the absence of 
an adequate heat sink to absorb the latent 
heat of condensation- It was suggested by 
the authors of reference 3 that the addi- 
tion to the flow of a noncondensable cool- 
ant gas should significantly increase the 
particle formation efficiency. It was 
therefore the purpose of this analysis to 
determine if the application of this tech- 
nique would result in particle formation 
efficiencies of 50$ or greater for a col- 
loid thrustor at realistic operating 
points. 

Over the past few years, a number of 
investigations have been conducted to 
evaluate the particle formation efficiency 
of a mixed flow condensation process 
(refs. 4 to 6). These approaches have 


been quite useful in developing a qualita- 
tive understanding of the condensation 
problem, but they have been of little help 
in the prediction of quantitative limits 
for particle formation efficiency as ap- 
plied to a colloid thrustor. 

Many of the analytical difficulties 
of the aforementioned approaches are 
avoided by basing the mathematical descrip- 
tion on an overall energy balance instead 
of details of the condensation process. 

Such a balance is obtained from the energy 
exchanges between the vapor, coolant gas 
and condensed particles. In this approach 
the latent heat rejected by vapor conden- 
sation is absorbed by several energy sinks 
in the system. These sinks for realistic 
thrustors are: 

(1) the thermal energy of the coolant 

gas, 

(2) the volume energy" of the con- 
densed particles, 

(3) the directed translational ' energy 
of the stream, 

(4) the thermal energy of the uncon- 
densed vapor in the stream. 

To further simplify the analysis, the con- 
densable vapor is assumed to remain in a 
saturated state during the condensation 
process. This condition does not appear 
to be too far from physical reality for a 
number of propellants in light of the ex- 
perimental data of references 5 and 7. 

PHYSI CAL MODEL 

A schematic diagram of the colloid 
thrustor analyzed herein is shown in fig- 
ure 1. A mixture of condensable vapor 
and noncondensable coolant gas is injected 
into a stagnation chamber where it is 
heated before being expanded in an adia- 
batic axisymmetric supersonic nozzle. 

During expansion, the vapor condenses to 
form small particles which, for a colloid 
thrustor, should be about 10^ atomic mass 
units per particle (ref. 8). The parti- 
cles then pass into the charging chamber. 
(The charging may take place by electron 
attachment as indicated in figure 1.) 

Then the negatively charged particles are 
electrostatically accelerated from a neg- 
ative thrustor potential to a potential 
slightly higher than that of free space. 
Finally the beam is neutralized and ap- 
proaches space potential. The nonconden- 
sable gas and any uncondensed propellant 
are lost to space in this concept . 

The overall thrustor efficiency is 
the product of the efficiencies associated 
with each of these processes (ref. 8). 
Associated with the nozzle condensation 
process is the particle formation effi- 
ciency, which is defined as the ratio of 
the condensed particle mass flow rate to 
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the total mass flow rate. This study is 
concerned only with the particle formation 
efficiency and the constraints placed upon 
it by the performance requirements of cer- 
tain components in the thrustor system. 

A schematic description of the phys- 
ical model used for the particle formation 
process is shown in figure 2. A mixture 
of vapor and coolant gas is heated in the 
stagnation chamber to some superheated 
state relative to the condensable vapor. 

The mixture is then expanded isentropic- 
ally until a saturated state is reached 
for the condensable vapor. In this anal- 
ysis, both the location of the saturation 
point, station 1, and its corresponding 
value of superheat required in the stagna- 
tion chamber, station 0, are fixed by 
specifying the value of the free stream 
Mach number at the saturation point ,M S . 

The nozzle expansion process from 
station 0 to station 1 is called region I 
of the particle formation process. This 
part of the expansion is shown as a solid 
line in both the temperature and pressure 
graphs. The form of the equations describ- 
ing the flow process in region I are the 
same as those for the isentropic expansion 
of a single gaseous component (ref. 9). 

The equation for the specific heat ratio 
and molecular weight of the mixture, how- 
ever, can be written in a form similar to 
that of reference 9 as. 
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All symbols are defined in appendix A. 

Expansion of the stream beyond sta- 
tion 1 results in condensation of the vapor 
phase. The maximum amount of condensation 
for a continuous condensation process oc- 
curs for an equilibrium saturated expan- 
sion of the condensable vapor. This con- 
dition is used to determine the stream 
temperature-pressure relation. The total 
heat capacity of the stream then governs 
the amount of propellant that will condense 
at any point in the nozzle. The experimen- 
tal data of references 5 and 7 indicate 
that it is possible for the vapor phase in 
a mixed flow nozzle condensation process to 
approach the conditions of a saturated ex- 
pansion. It is also indicated in refer- 
ence 7 that the amount of material that 
actually condensed was slightly less than 
the saturated value. 

The saturated expansion condensation 
process is considered to continue until the 
vapor phase becomes depleted or, until the 
flow leaves the nozzle. For the purpose of 
this analysis, vapor phase depletion is 
considered to have taken place when the 
mass flow rate of the remaining "uncondensed" 


vapor is equal to only 1 % of the total 
mass flow rate. Sample calculations indi- 
cate that an order of magnitude variation 
of this limit would not significantly af- 
fect the overall results of this analysis. 

The saturated expansion from station 1 
to station 2 is called region II of the 
particle formation process. The stream 
static temperature and pressure profiles 
are shown as dashed lines in this region 
in figure 2. Shown also for comparison, 
is a continuation of the initial isentropic 
expansion process from region I. 

Once the cutoff condition for the con- 
densation process is reached at station 2, 
it is assumed that any further expansion 
takes place isentropically until the noz- 
zle exit pressure is reached. As will be 
discussed later, the choice of the nozzle 
exit pressure must be compatible with the 
values required by the thrustor system 
operating restrictions. 

The portion of the nozzle expansion 
process from station 2 to station 3 is 
called region III. In region III, the 
isentropic expansion lowers the stream 
pressure and cools the condensed particles. 
The calculation is terminated wherever . 
stream pressure equals the preset value of 
nozzle exit pressure. It is assumed that 
no further condensation or stream expansion 
takes place beyond the nozzle exit. If a 
free expansion is considered beyond the 
nozzle exit, the condensed particle beam 
will probably become depleted of coolant 
gas and residual vapor in a short distance 
downstream of the nozzle exit due to the 
relatively high random velocity and low 
density of these species. Expansion of the 
condensed particle beam itself does not 
appear too feasible since the ratio of the 
random to directed velocity of the con- 
densed particles is relatively low. To 
collect and recirculate both the conden- 
sable vapor and coolant gas would be quite 
a formidable task for the operating con- 
ditions of interest. Therefore for the 
purposes of this analysis, all the residual 
vapor and coolant gas are rejected to 
space. Only the particle beam enters the 
charging chamber. 

To keep the analysis physically real- 
istic but analytically tractable, several 
simplifying assumptions are made. The flow 
is assumed to be one dimensional. Also it 
is assumed that a negligible amount of heat 
is transferred across the boundary layer. 
This physical model appears to be in agree- 
ment with the data, presented in refer- 
ence 10, in which stagnation and nozzle 
exit pressures similar to those of interest 
in this analysis (10 to 0.01 torr) were 
considered. In addition, it was shown in 
reference 1 that for operating conditions 
of interest to a condensation-type colloid 
thrustor, drop cooling by thermal radiation 
is negligible. Therefore, the entire noz- 
zle expansion process can be considered to 
take place adiabatically . Finally it is 
assumed that there is negligible condensa- 
tion on the nozzle walls so that continuity 
of mass flow is maintained in the stream 
itself. 
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ANALYSIS 


( 4 ) 


In this section, expressions for the 
stream properties in region II are ob- 
tained in terms of the initial conditions 
in the stagnation chamber and the propel- 
lant properties. 

Stream Temperature-Pressure Relation 

As discussed in the physical model 
section, the stream properties are control- 
led by the fact that the condensable vapor 
expands in a saturated state. The most 
common method of expressing the temperature 
vapor pressure relation, is by the 
Clausius-Clayperon equation 

dp v _ Tj- £v r 1 \ 

dT “uf v T \ l) 


^R_T 

P « (1 - g - f + fa) -y- 


Using equation (4) to eliminate \|r, one 
obtains equation (3) in the form 



Energy Equation 

Carrying out an energy balance between 
the stagnation chamber and any point in 
region II results in. 
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An expression for p y in terms of the 
stream properties can be obtained by using 
Dalton* s law of partial pressures with the 
equation of state and the continuity equa- 
tion. Three simplifying assumptions are 
made. First, the partial pressure exerted 
by the condensed particle gas is assumed to 
be negligible. Secondly, both vapor and 
coolant gases are treated as perfect gases. 
Since the pressure range of interest is 
relatively low, the compressibility factor 
for both condensable vapor and coolant gas 
is very close to unity. Thirdly, it is as- 
sumed that all species have the same di- 
rected velocity and temperature. Trial 
calculations Indicate that over most of the 
operating range of interest, the number of 
gas-particle collisions are sufficiently 
large so that this last assumption appears 
valid. However when the pressure ap- 
proaches the lower end of the range of 
interest (0.1 to 0.01 torr) this last as- 
sumption becomes marginal. The desired 
relation can now be written as. 
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Substituting equation (2) and the equation 
of state for the condensable vapor into the 
Clausius-Clayperon equation and differ- 
entiating yields. 



Momentum Equation 

The differential form of the momentum 
equation can be written as, 


f = - U du (3) 

Once again, if Dalton’s Law is used with 
the equation of state for each specie and 
the continuity equation, it is possible to 
obtain a relation between the stream pres- 
sure and density. 


Since, for the particle mass of interest 
(lO 5 amu), the random translational rota- 
tional and potential energies are negli- 
gible, the specific enthalpy of the con- 
densed phase can be written as 


C p T 

h C jf 
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Substituting this expression into equation 
(6) and differentiating yields, 

r— — « 
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Particle Formation Efficiency 

Equating the right sides of equations 
(2a) and (5), substituting the expression 
for U dU from equation (7) into the re- 
sulting equation ana then rearranging, 



Integrating this equation from the onset of 
condensation at the saturation point (sta- 
tion 1 in figure 2), to any point in 
region II yields. 



As a check on equation (9) it may be noted 
that when f is zero, the resulting equa- 
tion gives the condensation efficiency for 
a saturated expansion of a single vapor 
species (ref. 5). 

In order to obtain an expression for 
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the particle formation efficiency based on 
stagnation chamber conditions, a few modi- 
fications are required. The latent heat 
can be eliminated by using Troutons rule 
(ref. 11). 


92(3T b 
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the corrective term 3 is included to ac- 
count for the fact that some materials 
deviate slightly from Troutons rule. For 
most materials, 3 falls in a range between 
0.8 and 1.2. 

Using equation (10), the Clausius- 
Clayperon equation, the pressure and tem- 
perature relations for an isentropic expan- 
sion in region I and the perfect gas rela- 
tions, equation (9) can be rearranged to 
yield the desired expression for particle 
formation efficiency: 


condensation region, values of g and 
T/Tq must be found that satisfy both equa- 
tions (ll) and (13) . In the range of in- 
terest, 0 < g < (l - f), there is only one 
set of g and T/Tq that simultaneously 
satisfy both equations (ll) and (13). 

Since both these equations are transcenden- 
tal, they are not amendable to direct solu- 
tion by a simple hand calculation. The 
Newton-Raphson technique was used here to 
obtain a computer solution. 

LIMITATIONS 


Nozzle Flow Restrictions 

In order for the condensation process 
to take place, it Is necessary that a num- 
ber of flow restrictions be met. First, 
acting within the boundary layer are vis- 
cous forces that tend to dissipate the di- 
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Steam Pressure 

When Trouton's rule (eq. (10)) is used 
to express the latent heat in the Clausius 
Clayperon equation (eq. (l)), the result- 
ing integrated equation is 

f li.liPTv, \ 

p v = exp-j ? a + 11.513 + ll.lipf 

( 12 ) 

Normalizing equation (l2) to the condensa- 
ble vapor pressure at station 1, then using 
equation ( 2 ) and finally using the equa- 
tions for an isentropic expansion to ex- 
press the change in stream properties be- 
tween the stagnation chamber and station 1 
yields 


X ( 1 f ' ^ i r- + f a )exp ||ll . 513 + 


11.513+ 11.113 
-r/r-i 



(13) 


It should be noted at this point, that 
to obtain a value for condensation effi- 
ciency at any given stream pressure In the 


rected fluid energy so as to raise the 
fluid temperature. As a result, if there 
are any condensed particles in the bound- 
ary layer, they may be reevaporated by this 
viscous heat generation (refs. 5 and 12). 

An exact calculation of this effect is 
mathematically difficult, if not impossi- 
ble, to describe. Therefore an approximate 
approach is used. 

As a check to see whether or not the 
particles can be condensed in the boundary 
layer under operating conditions of inter- 
est for this analysis, the characteristic 
time for evaporating the particles can be 
compared with the nozzle residence times. 

As shown in appendix B, this characteristic 
time is on the order of 10 microseconds. 
When compared with the nozzle residence 
times which are on the order of 1 milli- 
second, it may be reasonably concluded that 
the condensed particles can be evaporated 
in the boundary layer. 

From the discussion above it can he 
seen that it is desirable to have the 
boundary layer occupy a small portion of 
the total nozzle cross sectional area. As 
a very optimistic estimate of where to cut 
off the nozzle condensation process due to 
boundary layer restrictions, that point in 
the nozzle where the boundary layer com- 
pletely envelops the flow was chosen (i.e., 
6/R - 1.0) . 

To locate the point in an axisymmetric 
nozzle where 5/R = 1, the empirical equa- 
tion 
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6 = Mh ! (14) 

1 yR? 

from reference 13 was used. In order to 
evaluate K, experimental values of 6 /l 
were chosen from reference 13, for Reynolds 
numbers (200 < Re < 7000), Mach numbers 
(3.5 < M < 6.0) and nozzle pressures 
(l0“2 torr < p < 10"^ torr) over the range 
expected in mixed flow colloid thrustors. 

A value of K * 0.7 was obtained. 

In terms of nozzle radius, equation 
(14) can be rewritten as 


Lan 1 / 2 1 


5-°- 7 .(f) - 1 . 

M 1 - 6 




Since A/A* , M and Re can all be related 
to the nozzle exit pressure, equation (15) 
expresses the relation between the nor- 
malized boundary layer thickness and the 
nozzle exit pressure for a fixed set of 
operating conditions and propellant prop- 
erties . 

As indicated in reference 14 if the 
flow is in either the slip or free molec- 
ular flow regime there will be a negli- 
gible amount of condensation due to the 
lack of sufficient expansive cooling. If 
the Knudsen number is used as a criteria 
for the establishment of continuum flow, 
it is found that continuum flow will take 
place when the Knudsen number is about 10' 2 
or less. Therefore, for the purposes of 
this analysis, the condensation process 
was considered to be complete when the 
stream state points are such that the 
Knudsen number is on the order of 1CT 2 . 

For the mixed-flow condensation process, 
this was found to be the case when the 
nozzle exit pressure reached a lower value 
of about 10-1 torr. 


Space Charge Restriction 

Since the colloid thrustor is an elec- 
trostatic device, the maximum flow rate is 
limited by space-charge effects. There- 
fore, all operating parameters such as 
stagnation pressure and nozzle exit pres- 
sure must be chosen to be compatible with 
space charge restrictions. To establish a 
criteria for the selection of these param- 
eters, the choked flow nozzle equation 
(ref. 9) Is combined with Child's Law for 
simple diode acceleration (ref. 15), and 
It is assumed that all the condensable 
vapor is condensed and all particles are 
uniformly charged. These assumptions re- 
sult in the following equation, 

• * «!/2 « pT/ r+i/r-i 

m sc = 8 € qRq e A 2 / r+ l\ 
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As discussed in reference 14, to comply 
with thrustor operating conditions, the 
maximum value of E A is about 4X10 6 volts 
per meter while the minimum value of spe- 
cific Impulse of interest is 2000 seconds. 
The choice of an applied field strength of 
4X10 6 volts/meter was based upon high volt- 
age breakdown data between flat parallel 
electrodes in vacuum. It is felt that this 
is a most optomistic choice, since in the 
presence of charged particles, residual 
vapor and coolant gas in the colloid thrus- 
tor accelerator, the maximum allowable ap- 
plied electric field will be somewhat 
lower. In addition, from practical consid- 
erations, the maximum value of stagnation 
temperature is limited to 1500° K. In- 
creasing this temperature as high as 2000° K 
makes little difference in the qualitative 
results. In order to meet space charge re-' 
strictions in an accelerator system, 
rfi sc > rhf . Therefore, the thrustor oper- 
ating conditions and propellant properties 
must be chosen so that the right side of 
equation (16) is always equal to or greater 
than one. 


RESULTS AND DISCUSSION 

As can be seen in equations (ll) and 
(l3 ) , there are an extremely large number 
of propellant properties and nozzle oper- 
ating conditions that can be varied to 
evaluate the particle formation efficiency. 
A number of these parameters, however, can 
be preselected based upon qualitative 
physical reasoning. Hydrogen was chosen as 
a coolant gas because it has the largest 
heat capacity. The large heat capacity is 
desirable because the amount of material 
that will condense Is limited by the heat 
capacity of the stream. The coolant gas 
can represent as much as 75 $ of the stream 
heat capacity in cases of interest. 

The choice of the value for Mach num- 
ber at the saturation point for the con- 
densable vapor, M s , is based upon two con- 
ditions. It is desirable to choose M s 
large enough so that the effects of thermal 
choking can be avoided. (As a result of 
thermal choking, all the material pre- 
viously condensed can evaporate (ref. 5). 
Sample calculations for a mixed flow sat- 
urated expansion indicate that when 
M s > 1.1 thermal choking can be effec- 
tively avoided. If M s is increased be- 
yond a value of 1.1 for fixed stagnation 
and nozzle exit pressure, the particle for- 
mation efficiency will drop since the ac- 
tual condensation process will occur over a 
lower expansion ratio. Therefore a value 
of M s * 1.1 is used as a limit. 
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Plotted in figures 3(a) and (b) is the 
particle formation efficiency as a func- 
tion of the fractional flow rate of cool- 
ant gas. Treated as a parameter in fig- 
ure 3, is stagnation pressure. As the 
value of the stagnation pressure increases 
from 1 to 100 torr, the maximum value of 
g increases from about 54 to 87$ while 
the optimum value of coolant gas flow rate 
drops from about 42 to 11$. It should be 
noted that these curves were drawn for 
fixed values of e =* 0.8, a = 30, 1.67, 

and P e = 0.1 torr. It was found that 
increasing € as high as 1.2 and dropping 
y as low as 1.3 tended to decrease the 
values of g for all cases to a slight 
degree. Varying a and P e had much 
more of a pronounced effect on the particle 
formation efficiency. Shown In figure 
3(b) is the effect of decreasing the noz- 
zle exit pressure to 0.01 torr. For -this 
case it can be seen that the maximum value 
of particle formation efficiency for 
P 0 = 100 torr is now about 90$ while that 
for 1 torr is about 78$. Thus, as might 
be expected, the particle formation effi- 
ciency increases with a decrease in nozzle 
exit pressure (all other parameters held 
fixed) due to the increased stream expan- 
sion. Although not shown in figure 3, in- 
creasing a tends to increase the values 
of particle formation efficiency. For ex- 
ample, for a stagnation pressure of 1 torr, 
a nozzle exit pressure of 0.01 torr and 
for values of other parameters as indi- 
cated in figure 3(b), the maximum values 
of particle formation efficiency are 
about 90$ and 58$ for values of a of 
100 and 10, respectively. Thus, if p 0 
and a are chosen high enough and p e 
chosen low enough, any desired value of 
particle formation efficiency should be 
attainable with no other constraints 
placed upon the system. 

For the case of a mixed flow colloid 
thrustor, however, there are a number of 
restrictions that must be met, as pointed 
out previously. Shown in figure 4 is a 
plot of 6/R as a function of nozzle exit 
pressure. Once again p 0 is treated as a 
parametric variable. Furthermore, this 
curve was drawn for the fixed values 
f = 0.3, e = 0.8, a = 30. It can be seen 
from figure 4 as the stagnation pressure 
increases from 1 to 100 torr, the minimum 
allowable nozzle exit pressure increases 
from about 0.1 torr to 0.6 torr. Stagna- 
tion pressures below 1 torr are of no 
interest since the nozzle flow Is no longer 
in the continuum regime. Experimental 
evidence (ref. 14) indicates that there Is 
a negligible amount of condensation when 
the flow is either in the transition or 
free molecular flow regimes. For those 
cases where the particle formation effi- 
ciency is greater than 50$, there Is no 
significant drop in the minimum allowable 
nozzle exit pressure for variations in 
f, and a around the values used in ob- 
taining figure 4. Therefore, it may be 
concluded that, due to the boundary layer 
restriction, the minimum allowable nozzle 
exit pressure Is on the order of 0.1 torr. 


The effect of the space charge re- 
striction on nozzle exit pressure is shown 
In figure 5. This figure is plotted for 
the fixed values f = 0.3, e = 0.8, and 
a = 30. From this figure it can be seen 
that for stagnation pressures between 1 and 
10 torr, the maximum allowable nozzle exit 
pressure is on the order of 0.01 torr. In- 
creasing the stagnation pressure beyond 
10 torr, drops the maximum allowable exit 
pressure well into the 0.001 torr range. 

To match the boundary layer restriction the 
highest possible nozzle exit pressure is, 
of course, desired. If the value a = 30 
is held constant and both e and f varied 
over wide ranges, the position of the 
curves in figure 5 do not significantly 
change. When a is increased above 30 the 
curves tend to drop to lower values, which 
Is in the wrong direction as previously 
pointed out. As a is decreased below 30,. 
the curves tend to rise, however if a falls 
much below 10, the particle formation effi- 
ciency is less than 50$. In general, it 
may be concluded that when 10 < a < 50 
and 1 < p 0 < 10, the maximum allowable 
nozzle exit pressure due to space charge 
restrictions is on the order of 0.01 torr 
and the particle formation efficiency is 
limited to about 80$. 

Another interesting result that can be 
obtained from the space charge restriction 
is that of the choice of a suitable cool- 
ant gas. As was stated previously, hydro- 
gen is desired as a coolant gas due to Its 
relatively high heat capacity. Also, equa- 
tion ( 16 ) shows that the ratio m so /mf 
varies inversely as the square root of the 
molecular weight of the coolant gas. In- 
creasing the weight much beyond a factor of 
2 would tend to drastically drop the curves 
shown in figure 5. As a result, the maxi- 
mum allowable nozzle exit pressure due to 
space charge restrictions would shift to 
much lower values than lO"^ torr. With hy- 
drogen as a propellant, the minimum allow- 
able nozzle exit pressure due to boundary 
layer restrictions is about 0.10 torr. 
Therefore by going to coolant gases with 
molecular weights much higher than that of 
hydrogen, an already difficult nozzle exit 
matching problem would become even more 
severe . 

Finally, shown in figure 6 is the par- 
ticle formation efficiency using the phys- 
ical properties of three possible propel- 
lants representing about an order of mag- 
nitude variation in molecular weight. In 
all cases the particle formation efficiency 
Is plotted as a function of the fractional 
flow of hydrogen coolant gas. Shown in 
figure 6(a) is a lightweight propellant 
(lithium), in figure 6(b) is a medium 
weight propellant (sodium), and in figure 
6(c) a heavv propellant (zinc). Comparing 
figures 6(a) through (c), as the propellant 
molecular weight increases, the particle 
formation efficiency also increases. 

In all cases these three specific pro- 
pellants examined encountered the same 
problems as the other general propellants 
previously examined in figures 3, 4, and 5, 
namely, it was not possible to match the 
minimum allowable nozzle exit pressure due 
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to boundary layer restrictions (0.1 torr) 
with the maximum allowable pressure due to 
space charge restrictions (0.01 torr). 

CONCLUDING REMARKS 

Over a wide range of operating con- 
ditions, propellants, and coolant gases, 
it was found that particle formation ef- 
ficiency was virtually unlimited. How- 
ever for cases applicable to colloid 
thrustors where the particle formation ef- 
ficiencies were greater than 50$, the min- 
imum allowable nozzle exit pressure due to 
boundary layer and non-continuum flow ef- 
fects were about 0.1 torr, while at the 
same time the maximum allowable nozzle 
exit pressure due to space charge effects 
was about 0.01 torr. Thus unique tech- 
niques to match the boundary layer and con- 
tinuum flow limits with the space charge 
restriction appear necessary. Expanding 
the colloid beam beyond the nozzle exit 
prior to entering the accelerator region 
appears to be impractical. 

APPENDIX A 

A area, m 2 

Cp specific heat at constant pressure, 

joules/kg 

E a applied electrostatic field, v/m 

reevaporation energy per unit volume, 
joules/m 3 

if frictional energy per unit volume per 
unit time, joules/m 3 -sec 

f fractional mass flow of coolant gas 

g particle formation efficiency 


P ratio of actual latent heat to that 

predicted by Trouton’s rule 

y gaseous specific heat ratio 

5 boundary layer thickness, m 

permitivity of free space, 
Q.85X10~12 C 2 /nm 2 

T] viscosity coefficient, nsec/m 2 

9 nozzle expansion half angle, deg 

p density, kg/m 3 

t time, sec 

t total dispersed density of mixture, 

kg/m 3 

Subscripts 

C condensed phase 

f nozzle flow 

n noncondensable coolant gas 

o stagnation chamber 

s saturated state at station 1 

sc space charge 

v condensable vapor 

Superscript 
* nozzle throat 

APPENDIX B 

The characteristic reevaporation time 
is defined as 



g 0 acceleration due to gravity on earth’s 
surface, m/sec 2 

h specific enthalpy, joules/m 3 

I S p specific Impulse, sec 

K constant 

L latent heat of vaporization, joules/kg 

l axial length of nozzle downstream of 
throat, m 

M Mach number 

JH molecular weight, kg/kg-mole 
m mass flow rate, kg/sec 

p pressure, n/m 2 

H radial distance to nozzle wall from 
nozzle axis, m 

He Reynolds number based on axial length 
R 0 gas constant, joules/mole-°K 

r radial distance from nozzle axis, rn 

T temperature, °K 

U directed velocity, m/sec 

U a directed velocity along nozzle axis, 

m/sec 

x distance, m 


where is the energy per unit volume 

required to completely reevaporate the 
condensed particles and /f is the energy 
dissipated per unit time per unit volume 
by the frictional forces. For steady 
state flow of a newtonian fluid in an 
adiabatic nozzle. 





2 
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Hi 
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(Bl) 


This equation, which was obtained by ma- 
nipulation of the Navier-Stokes, continuity 
and energy equations, is similar In form 
to that given in reference 16. In order 
to evaluate equation (Bl), it is necessary 
to know the spatial fluid velocity distri- 
bution . 

Since, over the operating conditions 
of interest, the Reynolds number for the 
mixed flow colloid thrustor nozzle is so 
low (5000 < Re), the flow can be treated 
as laminar. If the magnitude of the fluid 
velocity is assumed to vary parabolically 
with radial distance from the nozzle axis. 


if 


2 2 


4UU a r 


(B2) 


cl ratio of condensable vapor to coolant 
gas molecular weight 


In general, for laminar flow through an 
axisymmetric duct, a parabolic velocity 
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distribution occurs when the boundary layer 
completely envelops the nozzle. For the 
purposes of this analysis, it is assumed 
that no appreciable change in the particle 
formation process occurred due to fric- 
tional energy dissipation until this con- 
dition was reached. (This assumption 
yields an optimistic value of condensation 
efficiency . ) 

Combining the expression for the vol- 
umetric energy required to completely 
evaporate the condensed particles at any 
given point in the nozzle, 

» ^gL (B3) 


with equation (B2), 
relation 


x - 


results in the desired 


4TTU§r 2 

tgLR 4 


(B4) 


When equation (20) is evaluated over 
the range of propellant properties and 
operating conditions of interest, the 
characteristic times to completely evap- 
orate the condensed particles is on the 
order of 10 microseconds. The nozzle 
sizes chosen for this calculation corre- 
sponded to thrustor power levels as high as 
10 4 watts and specific impulse levels as 
low as 2000 seconds. Increasing the power 
level to as high as 10° watts does not 
qualitatively change these results. 
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Figure 1. - Schematic of a mixed-flow condensation type colloid thrustor. 
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(a) Static stream temperature against 
axial distance. 



ib> Static stream pressure against axial 
distance. 



(c) Particle generator. 


Figure 2. - Schematic description of the particle formation 
process in a mixed-flow, condensation type colloid thrustor. 
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(b) Nozzle exit pressure, 0.01 torr. 

Figure 3. - Particle formation efficiency as a function of fractional 
flow rate of coolant gas. 



Figure 4. - Effect of boundary layer restriction on nozzle exit 
pressure. 
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Figure 5. - Effect of space charge restriction on nozzle exit 
pressure. 



(c) Zinc propellant. 

Figure 6. - Particle formation efficiency for a mixed flow colloid 
thrustor. 


